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The time— dependent, spherically symmetric, Wyman sector of the Unified Field Theory 
is shown to be equivalent to a self-gravitating scalar field with a positive-definite, repul- 
sive, self-interaction potential. A homothetic symmetry is imposed on the fundamental 
tensor, and the resulting autonomous system is numerically integrated. Near the criti- 
cal point (between the collapsing and non-collapsing spacetimes) the system displays an 
approximately periodic alternation between collapsing and dispersive epochs. 

1. Introduction 

The Unified Field Theory (UFT)0i has held varying degrees of interest over the past 
50 years, however few would argue that the the essential idea (that of describing 
more physics by the "geometric" part of the field equations) was unwarranted. 
One of the interesting features of such theories is the inability to identify a purely 
gravitational sector; gravitation is enfolded in a larger structure that was originally 
intended to describe both gravitational and electro-magnetic phenomena. This is, 
of course, the reason for the name "Unified Field theory" , and inevitably leads to 
a theories for which it is difficult even to formulate a strong equivalence principle. 

Recently there has been some interest in claimSj-jtiiat a reinterpretation of the 
UFT, the nonsymmetjip gravitational theory (NGT),LfQ-|ai»ids singularity formation 
in collapse scenarios. Q^H This is-jppativated by the studyQa^u of the static, spherically 
symmetric, Wyman solution.t30 Claims that NGT should possess some repul- 
sive force that will overcome the attractive force of Einstein gravity in strong field 
regimes apd thereby prevent collapse to a black hole, have also appeared in the 
literature.ll3 This drew counterarguments, taking the form of 'thin shell' resultgl3 
and the contention that for initial data that is 'ekise enough' to a GR configuration 
the formation of a black hole is to be expected.EJ 

In this work we add somewhat to the debate by mapping the Wyman sector 
field equations to those of General Relativity coupled to a self-interacting scalar 
field (referred to as the Einstein-Klein Gordon (E-KG) system) in Section 2, show- 
ing explicitly that the potential is indeed repulsive with respect to the centre of 
symmetry. This shows up as explicit coordinate dependence in the potential, al- 
lowed in this case by the lack or separate conservation laws for the matter and 
gravitational fields (reflecting the lack of a strong equivalence principle) . We there- 
fore expect to see some signs of a slowing of the collapse of a spherically-symmetric 
matter distribution, and indeed, we show that for the simplified spacetimes with a 
homothetic symmetry considered in Section 3 (mimicking the results on collapsing 
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self-similar and homothetic E-KG spacetimeJli which may be important f or ii tljL e 
near-critical scalar field collapse that has appeared in numerical investigationst£lllII) , 
we see just such an effect. In this way we hope to gain some insight into what to 
expect from ihc more general numerical collapse problem which will clarify this 
issue further 

Perturbative analyses have shomuthat the original UFT/NGT field equations 



suffer from bad global asymptotie The most straightforward alteration ot-tbe 

action that avoids this criticisrrO'tSt^ is unstable near any GR configuration,E3H 
and is therefore untenable as a physical theory. Nevertheless, the spherically sym- 
metric Wyman sector considered alone suffers from none of these flaws, and has 
remained essentially unchanged throughout these alterations. Therefore we expect 
that when an action for NGT is discovered that is free of these problems (one such 
proposal existsc^) , the Wyman sector discussed herein will describe at least part of 
the spherically symmetric sector of the theory. 

2. The Wyman Equations 

To begin, we will review the field equations for a spherically symmetric, self- 
gravitating, self-interacting scalar field in Bondi coordinates. This will allow a 
straightforward map of the Wyman equations to this E-KG system with a specific 
choice of interaction potential. Not only does this simplify the analysis of the field 
equations, but many results from E-KG may be ported directly to the Wyman 
system (the identification of a horizon, causal structure, mass function, etc.). 

2.1. Spherically Symmetric, Self Interacting, Scalar Field 

Generalising the development of Christodoulou0 (using the same partially null 
frame) to include the self-interacting scalar field, we consider the spherically sym- 
metric spacetime metric in Bondi coordinates: 

=e^''du <»du + e^^^idu <^dr + dr<^ du) - r^dO ® dO - s\n^{9)d(j) ® dcf), (I) 

where u is a null coordinate which is constant on the future light cone of all points 
on the central world line, and r has been chosen so that spherical shells on these 
surfaces have area Anr'^. Furthermore we introduce the partially null frame 

en ■■= e'^du - ^e^^dr, er.^e'^dr, eg:^r^^de, := (rsin(6'))"^90, (2) 

in which the nonzero metric components are g^^^ = gf^^ = —ggg = — g^^ = 1- We will 
also make use of the standard definitionsEl / e''^^, g :— e''"'"'^, the mass function 
TO |(1 — e^^^) = |(1 — f /g), and the derivative along the ingoing radial null 
lines D := du — \fdr. 

The standard interacting scalar field Lagrangian {Vltjj] > 0) 

C ■■= lg''''d^[tP]dM - V[tP] - e-(^+^)9,[V;]i?M - V[i;], (3a) 

results in the stress-energy tensor for the scalar field T^j/ = 9^ [t/j] i9,y [t/'] — g/iu^, 
which has frame components 

T„„ = c-2-(Z?[V'])^ Tu^e-^^{dr[4>]f, Tnl = T,n = V[lP], Tgg ^ T^^ ^ C. 

(3b) 
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The frame components of the Einstein tensor 

Gf, =Rfi = ^e-^'drW + X], (4a) 

Gil =Ree = -^-^^^'^drire--'] " e^+'), (4b) 

^Kn = ^r^-''drW + A] - A) Q^JA] , (4c) 

G^ee ^Ki = -e-(''+^)a.[a„[z. + A] - e'^-^drM] + le-^^dr[iy - A], (4d) 

complete the description of Einstein's equations (IGttG = c = 1): G^^^ ~ i^m" ■ 
From the II and nl field equations respectively, we find 

dr[Hg)]=lr{dM)^, dr[rf]^{l~yV[^])g, (5) 

the former leading to the integral solution (assuming that ■(/; ~ i or faster as r ^ oo) 



9 = exp 



{d,ms)]fds, (6a) 



indicating that 5 is a nondecreasing function of r on each surface of constant it, 
and since we have fixed 5 ^ 1 as r ^ 00 for asymptotic flatness: < 9 < 1. The 
integral for / is 

f = g + f, where f ■= - J 5(5) yV^iV-K^) rfs. (6b) 

Here we have introduced the radial average hy g :— i Jj g{s) ds, and since g is 
a nondecreasing function of r we find that g < g (using g{u,r) < 9(11, rmux) for 
r G [0,rinax] in the integral for g). We have also assumed that / ^ 1 as r ^ 00, 
and that ij; and the potential V[ip] are sufficiently well-behaved so that /' ^ as 
r 00. Note that the presence of the potential results in /' > 0, and therefore 
/ > allowing the possibility that f > g and therefore m < 0. 
Since ingoing, radial, null geodesies are determined from 

du[xiu;r)] = -y{u,xiu;r)), (7) 

with solution x(u;ro) ~ ro — ^ J^^ /(s, x(s; ^0)) ds, the local magnitude of / con- 
trols the 'rate' at which the geodesic is proceeding towards r = 0. Defining the 
displacement vector between two nearby geodesies as C{u;ro) := drg[x{u;rQ)], we 
derive the result that ({u;ro) — exp[— i J^^ dr[f][s,x{s;ro)'^'^, and we see that the 
local focusing of the ingoing radial null geodesies is governed by 

dr[f]^l{9^f) = H9-9~n- (8) 

In the absence of the potential term we know that f ='g < g and so the geodesies are 
always focused towards r = (since then C < 1), however if f > g then defocusing 
may occur. 

It will also be useful to have an energy conservation equation; using 9m [m] — 
re"2Aa„[A] and (|), we have dr[m] = ^e-^^dr[y -t- A] -f \e~(''+^)r^V[^] and the 
nn equation may be seen to be an evolution equation for the mass function along 
ingoing radial null geodesies 

^gD[m] = ~UDm^-\fV[^]. (9) 
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Inserting into the Ggg component of the Einstein tensor yields 

G^ee = y-'drmDm - hn^] - irdr[V[^P]] - \rdr[^]V'[ij], (10a) 
which is equivalent to the remaining Einstein equation Ggg = \Tgg only if 

dr[<Pmi^]+dr[V[<P]\^Q, (10b) 

where the covariant wave operator on -0 is given by 

V^iV-] g^-^V^V.iV^] = lg'^D[dr[ri:]] - g-^dr[f]dM. (11) 
The equation of motion for ip derived from the action via the variational principle 

V^[i^]+5^,[V[^]]=0, (12) 



combined with (lOt) requires that dr [V[i)\\ = 5^ \V[ip]] drW\, which is true provided 
there is no explicit dependence on the radial coordinate in the potential. We will 
have cause to re-examine this in the following section. 

2.2. The Wyman Equations 

We will avoid giving the details of the computation of the field equations (the mote 
general reduction of the system to (1 + l)-dimensions wilL appear elsewherecS) 
instead quoting results derivecj— ^^sjng an existing method.Ej The field equations 
may be derived from the actio: 

S = 



(13a) 

where 

The Wyman equations are those that result when the spherically symmetric ansatz 
(the functions 7, (5, a, /?, and / are functions of i and r only) 

g =^dt iS) dt + S (dt dr + dr dt) — adr ® dr 

(14) 

- I3{d9(g)d0 + siii^{0)d(l)(g)d(l)) + fsm{e){d0(g)d(j)-d(l>(g)de), 

for the fundamental tensor is inserted iato-thc. general field equations, and have 
appeared various places in the literatureJiO'Ej'Ej'Ej 

Since the t-r sector of (^^ is symmetric, we can once more make use of Bondi 
coordinates, and furthermore parameterise the antisymmetric part of the Wyman 
tensor using a dimensionless angle 4' as 

g =e^''du ®du + e^^^^du ® dr + dr ® du) 

- r^ cos{tp){dd (g)dO + sm'^{9)d(j) (g) d(p) - sin(?/') sm{6){d9 dcj) - dcf) ® dO), 

(15a) 
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which becomes 



cos 



3in(V')(6i^0 6''^-6''^i 



(15b) 



This parameterisation is 
(i.e., is identi cal t o the volume element 



in the partially null frame introduced in Section 2.1. 
particularly convenient since ^/— g = 

of t he m etric (n])) and the components of the inver se of (15b) are identical to that 
of (15b) with -ip-^ —ip (the angular components of (15b) have the form of a rotation 
matrix) . _ 

The field equations (in the form = CO) may be found by straightforward 
computation, and are related to the results of the last section by 



0, 



where / and g are defined as in Section 2.1, and the potential is given by 



1 



3(V')) + iM'sin2(V). 



(16a) 
(16b) 
(16c) 
(16d) 
(16e) 

(17) 



For explicitness we give here the field equations in the form most useful for later 
considerations: 



dr[rf] = cos(V')5 - jiprfgsm^itp), 



(18a) 
(18b) 

2rD[dr[r^/j]] - dr[f]dr[^j] + 25sin(V') + (^r)^^ sin(V') cos(V') = 0. (18c) 

It is a rather laborious task to show that the Wyman solutio identically 
satisfies these field equations with u = 0. 

The first three equations of (|l6[) ^re clearly equivalent to the ll^nl and nn 
field equations of the E- KG s ystem (|) with the potential given by ([l7|), leading 
to (^a]), (|6b|) , and M). In (16d) we fina the field e qu ation for i/j as woulobe de rived 
from the variational principle with the potential (O). The final field equation ( |l6e ) 
is somewhat odd looking since it would appear to oe incompatible with both of ([17) 
and the field equation Ggg — ^Teg. However if we allow for the possibility that the 
potential depends explicitly on the radial coordinate, then the partial derivative of 
the potential splits up into two parts (the second term represents the derivative of 
the potential with respect to the explicit coordinate dependence only) 



(19) 



and instead of (10a) we have 



\rdAVM ~ \rdrm{V^[^P]+S4V[^P]]). (20) 



However by direct computation, we finds from ( |l7| ) that 



(21) 
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and therefore the field equations are consistent. 

As advertised, the potential (O) acts (for non-zero ^p) as a repulsive central 
force; indeed we expect that it wiirovercome the attractive gravitational potential 
at small enough radii. In the E-KG system, explicit coordinate dependence in the 
scalar potential is in conflict with the Bianchi identities; in flat spacetime leading to 
conservation laws of the formE3 9 • T = —dC. rA4fe^tress though that what is playing 
the role of the Bianchi identity in UFT/NGTBB is not the same as in GR. For the 
E-KG system we had (for spherical symmetry) a field equation for the scalar field 
derived from the variational principle that was equivalent to the conservation laws 
(and therefore Einstein's field equations) provided there was no explicit dependence 
on the coordinates. Here we have a (diffeoniorphisni invariant) system of equa- 
tions that includes the scalar field equation that altogether satisfies the UFT/NGT 
Bianchi identities, but does not satisfy the separate E-KG conservation laws. 

For this work the operationally important point is that the functions that pa- 
rameterise the symmetric part of the fundamental tensor are determined in exactly 
the same way from the scalar field as if it were a spherically-symmetric, E-KG field. 
Furthermore, the field equation for ^ is exactly that which would be determined 
from the variation of the scalar field action. Therefore we can treat the system as 
E-KG, ignoring the fact that the potential has explicit coordinate dependence. 

3. Self Similar Solutions 

Evidence for self-similarity in the numerical collapse of a scalar field was first noted 
by ChoptuikJl3 and sparked much interest in the literature. Since that time it 
seems accepted that it is a discrete self-similarity rather than therfehe continuous 
(or homothetic) self-similarity considered herein which is observed.Ea Nevertheless 
we feel that there is something to be learned in examining the homothetic solutions: 
not only are the resulting equations much simpler to integrate, but the behaviour 
which results is reflecting the repulsive nature of the potential (|^ and should be 
observable in more generic oeilapse scenarios. 

We now make the ansatzc3 that the fundamental tensor has a homothetic sym- 
metry, that is, there exists a vector field ( for which the fundamental tensor satisfies 

£ds\^2g. (22) 
It is straightforward to show that for the fundamental tensor ( p^ ) we have 

£c[X]=0, £cH=0, £cM=0. (23) 
where the similarity vector is 

C = udn+rdr, (24) 

and so all three functions must depend solely on the variable x := r/u. From the 
results of the last section we know that this is equivalent to considering a homothetic 
E-KG system with a nontrivial potential, restricted to the scale invariant case 
(k = 0) of Brady.B 

Indeed, normally if one considers a scalar field potential and requires spacetime 
to have a homothetic symmetry, one ends up with the requirement that i?^ [^[V']] + 
21/['(/)] = 0, which for potentials that are polynomial in the scalar field, requires 
that "0 = 0, and one is left with a trivjial spacetime. To see this we follow the 
argument in the appendix A of Bradyll3 and note that £Q[R]fii, = now implies 
that £l:[Vd^p]V4^|J] - gf,^V[:p]] = 0. Considering an angular component of this 
and using spherical symmetry and (|2^ ) yields the result. Using the potential (|l7| ) 
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and i^f ["0] = 0, we satisfy this requirement due to the exphcit dependence on the 
radial coordinate in the potential. Therefore the existence of self-similar Wyman 
spacetimes is not a trivial extension of the E-KG results, it is actually a rather 
surprising possibility that is wrapped up in the structure of the field equations. 
(Note that we are considering /i = throughout; the existence of this length scale 
would clearly contradict scale invariance. We expect however that if the homothetic 
symmetry is relevant to the general collapse problem, then the generie features 
presented herein should also appear for fi ^ on distance scales < 
Inserting the ansatz into the field equations (|l^) we find 

d^[\ng] = lx{dMy, (25a) 
d,[xf]=cos{^)g, (25b) 
2xix + i/)9f M + x''d,[f]dM " 2gsin(V') = 0. (25c) 

It will be important later to determine the large x asymptotic behaviour of asymp- 
totically flat solutions of (pq). Assuming that / — > /o and g —> go {i.e., spacetime is 
flat in t he asymptotic limit) we find that the asymptotically dominant contribution 
to ( ^5c| ) is 

x^d'^i^jj] + 2x^d^[^lj] - go sin(V') - 0. (26) 

If ■(/;—!■ V'o 7^ 0, we can approximate sin{ijj) ~ sin('0o) +cos(?/;o)V' where ?/> —>■ V-'o + "01 
and then taking |a;| —: gocos{ijjo)y and X := tp + tan('i/;o) we have 

y'dl[X]+2y'dy[X]TX-0, (27) 

where the upper sign is for a; > and the lower for x < 0. Introducing y — aj it 
is easy to see that the solutions to this are Besscl functions, and once we requires 
that i/jQ we have the asymptotic form of the solutions 

i^x>o -i^o + t&n{il;Q){zK^{z) - 1) + Czh{z), (28a) 
V'x<o -V-o - tan(V'o)(^7rzYi(z) + l) + CzJi(z), (28b) 

where C is an integration parameter. As it turns out, the case where 7p Q la 
equivalent to taking -00 = in these solutions. 

Returning now to the system (p5[), we introduce the alternate variables Y := f / g, 

Z := x/ f, the coordinate t = ln|x| and the additional variable P :— tp (an overdot 
indicating a derivative with respect to t), and find the autonomous system 

Z =Z{2-Y-^ cos(0)) , (29a) 

Y = cos(0) - r - iyp2 ^ ^|±jl (cos(V') - Y) , (29b) 

tp =P, (29c) 
P = - y"i(2Z+ l)"i[(2rZ + cos(0))P- 2sin(V')]. (29d) 

where in ( ^9b| ) we have made use of the energy conservation equation (^ which, 
assuming homothetic symmetry, may be written in the useful forms 

„ 4(1 -r-icos(i/;)) 2 , , 

p2 ^ _V = dtllnY]. 30 

2Z + 1 Z+1 ^ ^ ^ ' 
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Note that {— oo < t < 00} covers either w < or < u, however the sohitions 
are mapped into each other by t <-> — t. We will be considering the former case 
exclusively where spacetime is empty and flat for t < 1 (so that = 0, Y = 1 and 
Z — x/fi for some constant /;). By rescaling of the u coordinate fo may be chosen 
at will; however this merely shifts t to + 1, mapping each trajectory onto itself. 
At t — 1 the scalar field will be 'turned on' {i.e., matched to a self-similar solution 
through the initial data) and the system evolved as far as possible towards < = 00 
{u = 0). 

3.1. The Einstein-Klein Gordon Limit. 

The equivalent set of equations for a non-interacting, self-gravitating, spherically 
symmetric, scalar field are equivalent to the Vltp] = limit of (E9): 



Z =Z{2-Y-^), (31a) 
^=P, P = -Y-\2Z +1)-'^[{2YZ +1)P], (31c) 



and 



p2 ^ 4(1 Y 1) ^ (32) 

2Z + 1 Z + 1 *^ J' ' 

replaces ( |30|) . 

Since afree scalar field has a positive mass function by construction (see (6b) 

with V[ip] = 0) we know that Y <1, we furthermore confine ourselves to y > smce 

y = signals the presence of an apparent horizon. With this information in hand, 

1 

' 2 ■ 



the requirement that (32) be positive (and therefore P real) results in Z < — i. 
Ingoing radial null geodesies may be reparameterised in terms of x as x(u;r) — 
ux{u;x) which, inserted into (|^) and integrated from u = —1 (corresponding to xi) 
to u gives 

/■X 

ln|u| = -2/ dx' (fix) + 2x)~\ (33) 



and with f = x/Z and Z < we see that all of the region of interest is accessible, 
and either u — )■ while Z ^ 00 or u uq a.t the horizon Z = — -i. 



Taking the ratio of (31b) and (31a) we can write Y = Y[Z] to find 



dY _ 2{Z+1) Y{1 - Y) 
'dZ ~ Z{2Z+1) 2Y -I ' 



(34) 



which implies that trajectories may only pass through the lines Z = Y ~ 
or F = — 1 at either (— ^,0) or (— and must cross Y = \ pe rpendicularly 
unless Z = —1. The point (—1, i) is an unstable fixed point of ( ^jla]) and (|31b|) Jij 
(We refer to points in the Z-Y , and later Z-X, plane by their coordinates: (— O) 
refers to the point Z = — i, F = 0.) Additionally, from the second form of ( p^ ) 



y > implies that F < for -1 < Z < -i and F > for Z < -1, and (31a) 



implies that Z < for i < F < 1 and Z>OforO<F<i. Therefore trajectories 
that are matchable to a fiat spacetime (only C^) emanate from (--^,1) ending at 
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either (— ^,0) (a horizon), (—1,^) (asyniptoticaUy approaching the fixed point), or 
traveUing to infinite Z with asymptotic form 

Y^l + Yoe-\ Z^Zae\ (35) 

which approaches a fiat spacetime configuration. Nothing need be said of ip since it 
does not directly enter into the equations for Y and Z; however since P is of a fixed 
sign throughout these trajectories, if -0 = initially it cannot vanish again. For a 
free scalar field this is irrelevant since a constant ip solution has vanishing energy 
momentum tensor. This will clearly not be the case for the homothetic Wyman 
spacetimes. r-i 

It is useful that there exists an exact solution to this system,E3 which we will 
reproduce here. Integrating (^) results in 

Yil-Y) = --^^, ^ Y=—{Z±^Z^+A{2Z+l)), (36a) 

where A > is a parameter of integration chosen so that the critical solution 
corresponds to A = Ac := 1, sub-critical < ^ < 1 and super-critical A > 1. Note 



that the ± sign in (36a) refers to i < y < 1 and <Y < ^ respectively. Inserting 
this into (31a) and requiring that Z ~ att — gives 



1 



■ (a + v/^(l-e2*)+e2*) . (36b) 



Note that super-critical solutions reach ^,0) in a finite parameter time t = 
iln(A/(A-l)). 

In Figure ^ we compare this exact solution ( p6| ) (indicated bv diamonds) to that 
derived from numerically integrating the system of equations (|3l| ) (solid lines) , also 
indicating the exactly critical solution {A — 1). In order to avoi3~the poor behaviour 
of {2Z + 1)^^ in ( ^T|) wc begin the numerically generated trajectories near ^, l) 
by considering initial data from the exact solution (|3^) Taylor expanded to a small 
parameter time At 

Y^l-AAt, Zw-i(l + Ai), V~±2\/AAt, P^±2VA. (37) 

Note that the choice of sign is immaterial, merely replacing tp ^ —ijj and P ^ —P 
along the entire trajectory. The initial time At was chosen by AAt = 10^^ which, 
for the range of A considered throughout this work, was sufficient for considering 
qualitative features of these spacetimes. In generatiafcall trajectories, we made use 
of the numerical ode package of Maple V Release 4;L3E3 although we have performed 
no quantitative error estimates, the agreement of the numerical versus exact results 
in Figure |l|, qualitative agreement with the analytic conditions derived, and the 
robustness of the results using the different integration schemes supplied by Maple 
provide confidence in the results. The results presented were generated using the 
Livermorc Stiff ODE solver (Isode) routine with default options. 



3.2. The Wyman Equations 

As before we begin by noting that we must have F > in order to be considering a 
(part of) spacetime that is not in the interior of a black hole. In this case however, 
we cannot claim that Y < 1 since by (pb| ) we need not have a p ositive-definite mass 
function. However using the explicit form of the potential (|l7| ) (with /i = 0) we do 
find that f = ^ dr cos{ip)g, which allows the simple bound |/| < g, and we need 
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-1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 



Figure 1: E-KG trajectories. The lines are calculated by numerically integrating 
the equations (^) , whereas the diamonds are calculated via the exact solution ( |36| ) . 
The dotted line indicates the critical solution which asymptotically approaches the 
fixed point at (—1, 5)- 



only consider spacetimes with < F < 2. In spite of this, we shall see that all of 
the trajectories we consider will remain within the region < y < 1. 

It is useful to consider the new variable X :— Y j cos('(/'), in terms of which 
takes on a form which is similar to ( |3l| ) 

Z =Z{1- X-^), (38a) 

^ =^ItT^(^ " + XPtan(V^), (38b) 

ii =P, (38c) 

P = - X-\2Z + 1)-^ [{2XZ +l)P-2 tan(V')] , (38d) 

and 

P^='^'-^-'K-^0,[lnX]. (39) 

2Z+1 Z+1 ^ ^ ^ ' 

The energy conservation equation is now identical to (^), and since there is no 
physical bound on X, we find that trajectories must be in one of the following 
regions: 

X<0, Z > 

0<X <l, Z <-\ . (40) 



X>\, Z> -\ 
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-us ^.7 ^.6 ^.5 U_ii _i ^7 6 ^5 

z z 



(a) Wyman Trajectories in the Z-X (b) Wyman Trajectories in the Z-Y 

plane. plane. 

Figure 2: A selection of the Wyman trajectories showing sub-critical {A = 1 and 
A — 2.5), sub-critical approaching {A = 3.3) and crossing {A = 3.5) the line 
X = 0.5, super-critical on the second oscillation with the first oscillation not passing 
through X — 0.5 (A = 4.5) and passing through X = 0.5 {A — 5.2), and finally a 
'direct' super-critical trajectory {A = 6.5). 



From our initial conditions, we require that ^ — and Y = 1, implying that 
trajectories of interest must begin at (— For a suitably small interval in t 
(where the tan(V') term in ( |38b| ) has little effect) we will find approximately E-KG 
evolution, and therefore trajectories will flow into region ii. Moreover, since we 
will never find that ^|J ^ (n + ^)7r on any trajectory, we may also unambiguously 
identify X = as a horizon, and may therefore consider region ii exclusively where 
again we see from ( ^3[ ) that the entire region of interest is accessible to incoming 
radial null geodesies. The fixed point (—1, i) that existed for the E-KG system 

may no longer be considered a fixed point since we have P — ip — 4: ^ (note that 
P = A also for the E-KG system, however there is no potential to feed back into 
the {Z,Y) system; this Tnea.n s th at •;/; ^ ipQ + At along this 'fixed point'). 



Taking the ratio of ( 38b ) to (38a) we find 



dX 2{Z+1) X{l-X) X^ Ptan(V;) 
'dZ~Z{2Z+l) 2X-1 '^2X-1 Z ' ' ' 

As for the E-KG system, trajectories may not cross Z = or X = except at 
(— i,0) or (— i,f), however they may intersect X = 1. Trajectories that do so at 
Z 7^ — i must intersect with P = from (p9|), so if we consider initial points at 
X = 1, Z ^ — \ then we initially have ■0 = -P = 0. It is easy to see that the solution 
is given by -0 = 0, P = 0, X = 1, and Z — Zqc' which is the flat spacetime solution. 
From ( p^ we deduce that trajectories will pass through X = ^ perpendicularly 

unless P{Z -I- 1) = 2tan(-0), which implies that P = and therefore that X = 
and so are always perpendicular. 

We display a variety of trajectories emanating from (— ^, l) in Figure ^. All sub- 
critical trajectories eventually approach X — \ which, as we see from Figure ^, does 
not in general correspond to = 0; instead their form corresponds to that given 
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Figure 3: The scalar field tp for the trajectories displayed in Figure ^ as well as the 
two near-critical trajectories in Figure |4| note the logarithmic time scale. 



by (pSbl ). However since the asymptotic form does not correspond to a solution 
of (P^J, we cannot match onto a flat spacetime solution at u = (as one could for 
the Fr-KG case, taking ip = constant). It seems likely that the trajectory would 
continue on through u = 0, matching to the equivalent u > trajectory, effectively 
following the same paths as shown in Figure back to Y — 1, at which point we 
could match back onto Minkowski space. We also find that there are no traject orie s 
with ip ±7r/2 (those which approach ±i7r are repelled by the tan(?/') term in ( 38d ) 
since P « 0) passing through (—5,0), so we conclude that there is always a horizon 
and need not continue the trajectory into region i. 

As we see from Figure ^, there is a qualitative transition at A = « 3.34 
where trajectories pass arbitrarily close to the point (—1, 5)- Whereas trajectories 
with A < Ai display no tendency to collapse, those with A > Ai undergo at 
least one period in which the gravitational attraction overcomes the tendency of 
the potential (|l^) to disperse the scalar field. Trajectories with A < Ai will never 
cross X = i, those with Ai < A < Ac ^ 3.563 will cross X = ^ an even number 
of times, eventually propagating to Z —i- 00, and those with A > Ac will cross 
X = i an odd number of times, eventually ending at (— ^, O) . Near-critical solutions 
(A « Ac) shown in Figure |^ display an approximate periodicity, resulting from 
alternating periods in which either the gravitational attraction or the repulsive 
potential dominates. 



4. Discussion 

Understanding these results is rather straightforward. In addition to the presence 
of the self-interaction term, the repulsive nature of (O) is acting quite strongly to 
disperse the field. This is shown by the increase in tne amount of energy required 
to produce a black hole over that of the E-KG system; Ac — I for the E-KG system 
and Ac ~ 3.563 for the Wyman equations. We can also understand the oscillatory 
behaviour as the interplay of the attractive gravitational force and the repulsive 
potential. 
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Figure 4: Near-critical Wyman trajectories; the sub-critical trajectories {A — 3.4 
and A — 3.563) are shown as solid lines while the super-critical trajectories A — 
3.564 and A — 4.6 are shown as dashed lines. 



Trajectories corresponding to very small initial energies (A « 0) are very close 
to their E-KG relatives; the field does not get large enough to feel an appreciable 
repulsive potential. As we consider increasing initial energy (increasing A) the 
potential quickly manifests itself by 'flattening out' the trajectories in Figure ^(b) 
over those in Figure ^; trajectories that were super-critical for the E-KG system 
have been repulsed by the Potential (pTj ). Notice that the tendency of the potential 
to disperse the field also has the effect of decreasing the repulsive potential, so 
at some point the attractive gravitational potential will take precedence and the 
system will have a tendency to re-collapse. For trajectories with initial momenta 
that is well above the critical value, enough energy is directed towards r — 'fast 
enough' that a horizon forms before the repulsive potential has an opportunity 
to take effect (as for the E-KG system these are infinitely massive black holes, 
and the entire spacetime becomes trappedcj). Close to the critical trajectories, 
the system oscillates between attractive and repulsive phases many times before 
'choosing' its fate. In a more general numerical simulation,E3 this behaviour should 
be straightforward to identify. 

So far we haws not discussed the ambiguity in the identification of the spacetime 
metriG4ij-4JFT,Lj an issue that is far from clear given the existence of multiple light 
cones.eifEj For this system, the ambiguity vanishes for the following reasons: For 
purely spherically symmetric fields the light cone is determined solely in terms of 
the metric in the t-r sector, and, as usual, we only need to know this metric up to 
an arbitrary conformal factor. The ability to reduce the system to an E-KG system 
tells us that the (radial) causal structure determined from the t~r components of the 
fundamental tensor therefore has physical meaaing. (It is also not difficult to show 
that the metrics determined by Maurer-TisonEa are all conformally equivalent to a 
metric with the same t~r components discussed here.) This unambiguously provides 
the causal boundary for Wyman sector field propagation, and results from GR based 
on radial null geodesies (such as the existence of a horizon) are therefore relevant 
for this problem. 
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